Higgs-Boson Production Associated with a Single Bottom Quark in 

Supersymmetric QCD 



Hou Hong-Sheng 2 , Ma Wen-Gan 1 ' 2 , Wu Peng 2 , Wang Lei 2 and Zhang Ren- You 2 
X CAST (World Laboratory), P.O.Box 8730, Beijing, 100080, People's Republic of China 
2 Department of Modern Physics, University of Science and Technology of China (USTC), 
Hefei, Anhui 230027, People's Republic of China 



Abstract 

Due to the enhancement of the couplings between Higgs boson and bottom quarks in the minimal 
sypersymmetric standard model (MSSM), the cross section of the process pp(pp) — > h°b(h°b) + X 
at hadron colliders can be considerably enhanced. We investigated the production of Higgs boson 
associated with a single high-pr bottom quark via subprocess bg(bg) — > h a b(h°b) at hadron colliders 
including the next-to-leading order (NLO) QCD corrections in MSSM. We find that the NLO QCD 
correction in the MSSM reaches 50% ~ 70% at the LHC and 60% ~ 85% at the Fermilab Tevatron 
in our chosen parameter space. 



PACS: 14.80.Ly, 11.30.Pb, 12.60.Jv 



I. INTRODUCTION 



One of the most important tasks of future high-energy experiments is to search for scalar Higgs 
particles and investigate the symmetry breaking mechanism of the electroweak SU(2) x U(l). In the 
standard model (SM) pP, one doublet of complex scalar fields is needed to spontaneously break the 
symmetry, leading to a single neutral Higgs boson h . But there exists the problem of the quadratically 
divergent contributions to the corrections to the Higgs boson mass . This is the so-called naturalness 
problem of the SM. One of the good methods to solve this problem is to make supersymmetric 
(SUSY) extensions to the SM. Then the quadratic divergences of the Higgs mass can be cancelled by 
loop diagrams involving the supersymmetric partners of the SM particles exactly. The most attractive 
supersymmetric extension of the SM is the minimal supersymmetric standard model (MSSM) [21 El- 
In the MSSM, there are two Higgs doublets Hi and H2 to give masses to up- and down-type fermions. 
The Higgs sector consists of three neutral Higgs bosons, one CP-odd particle (A °), two CP-even 
particles (h° and H°), and a pair of charged Higgs bosons (H^). 

In the SM, the so-called Yukawa couplings describe the interactions between Higgs and fermions. 
Their coupling strengths are proportional to rrif/v, where v is the vacuum expectation value of the 
Higgs field. Since the bottom quark mass is approximately 5 GeV, the Yukawa coupling to bottom 
quarks is relatively weak. This leads to a small cross section for the associated production of the Higgs 
boson (h ) and bottom quarks at hadron colliders [3J0I]. However, in the MSSM, the Yukawa coupling 
between the Higgs boson (h . ) and bottom quarks can be considerably enhanced for large values of 
tan/? = V2/V1, where v% and V2 are the vacuum expectation values of the two Higgs boson fields H\ 
and H2, respectively. 

Because the high-p^ bottom quark can be tagged with reasonably high efficiency, the observation 
of a bottom quark with high px can reduce the backgrounds of the Higgs boson production. The 
leading-order subprocess for production of the Higgs boson associated with bottom quark is bg —* bh° 
U3- Recently, the investigations of the process pp(pp) — ► bg(bg) —* h°b(h°b) + X in the SM including 
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the NLO QCD corrections is presented in Ref. j7|. Their calculation shows the cross section of the 
subprocess of bg — ► h°b is an order of magnitude larger than those of gg(qq) — » at the Fermilab 
Tevatron and the CERN LHC. They find that the NLO QCD correction ranges from 50-60% at the 
Tevatron for = 100 — 200 GeV, and at the LHC the correction with m/j = 120 — 500 GeV ranges 
from 20 — 40% for pt > 15 GeV and 25 — 45% for pr > 30 GeV. They conclude that this production 
mechanism improves the prospects for the discovery of a Higgs boson with enhanced coupling to the 
bottom quarks. 

In this paper, we calculated the cross section for the associated production of the Higgs boson and 
a single high-py bottom quark via bg — > h°b(bg — ► h b) in the MSSM at the Tevatron and the LHC 
including the NLO QCD corrections. The structure of this paper is as follow: In Sec. II, we discuss 
the LO results of the subprocess bg — ► h°b. In Sec. Ill, we present the calculations of the NLO QCD 
corrections. In Sec. IV, the numerical results and conclusions are presented. Some lengthy analytical 
expressions are listed in Appendices. 

II. THE LEADING ORDER CROSS SECTION 

Since the cross sections for the subprocess bg — > h°b and its charge-conjugate subprocess bg —* h°b in 
the CP-conserved MSSM are same, we present only the calculation of the subprocess b{p\)g(p2) — ► 
hP(k%)b(kiC) here (where p\^ and £^4 represent the four-momentum of the incoming partons and the 
outgoing particles, respectively.). The subprocess bg — ► hPb can occur through both s-channel and 
t-channel as shown in Fig.l(A-B). So we divide the tree-level amplitude into two parts and denote it 
as 

M° = M^ s) +M^\ (2.1) 

where and represent the amplitudes arising from the s-channel diagram shown in Fig. 1(A) 
and the t-channel diagram shown in Fig. 1(B) at the tree-level, respectively. The explicit expressions 
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for the amplitudes Mq S ^ and Mq can be written as 

M (s) = gs(Mr) g r6(Mr) ^(fe 4 )^+^)7^ i (p 1 )6,( P2 )^ ) 

M (t) = S.0h)YM ^ . {kA)lv{i>1 _ h ) Uj { Pl )e v {p 2 )T^ (2.2) 

where s = (p\ +P2) 2 , t = (pi — k^) 2 and u = (pi — k^) 2 are the usual Mandelstam variables. g s (fi r ) is 
the running strong coupling strength and T a is the SU(3) color matrix. lfe(/U r ) is the Yukawa coupling 
between Higgs boson and bottom quarks. In MSSM, Yb(^ r ) is given as 

w \ .9wfn b (ii r ) sina 
2mw cosp 

where a is the mixing angle which leads to the physical Higgs eigenstates h° and H°. mi,(fJL r ) is the 
MS mass of the bottom quark. Throughout our evaluation we adopt the simplified Aivazis-Collins- 
Olness-Tung(ACOT) scheme (Hj, which will not make loss of accuracy in our calculation of subprocess 
bg — > h b. That means the bottom quark mass is neglected except in the Yukawa couplings during 
our calculation. 

Then the lowest order cross section for the subprocess bg — ► hPb in the MSSM is obtained by using 
the following formula: 

I f'tmax 

£0 (|j bg ^ h o b) = __ di ^| M °| 2 , (2.4) 



mi n 



where i max = and t m i n = m| — s. The summation is taken over the spins and colors of initial 
and final states, and the bar over the summation recalls averaging over the spins and colors of initial 
partons. 

III. NLO QCD CORRECTIONS 

The NLO QCD contributions to the associated production of the Higgs boson and a single bottom 
quark can be separated into two parts: the virtual corrections arising from loop diagrams and the real 
gluon emission corrections. 
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1. Virtual Corrections 



The virtual corrections in the MSSM to bg — > h°b consist of self-energy, vertex and box diagrams 
which are shown in Figs. 2-3. Fig.2 shows the one-loop diagrams of the SM-like QCD corrections from 
quarks and gluons, and Fig. 3 presents the one-loop diagrams of the SUSY QCD corrections from 
squarks and gluinos. There exist both ultraviolet (UV) and soft /collinear infrared(IR) singularities 
in the amplitude from the SM-like diagrams in Fig.2, and the amplitude part from SUSY QCD 
diagrams (Fig. 3) only contains UV singularities. In our calculation, we adopt the 't Hooft-Feynman 
gauge and all the divergences are regularized by using dimensional regularization method in d = 4 — 2e 
dimensions. 

In order to remove the UV divergences, we need to renormalize the wave functions of the external 
fields, the strong coupling and the h° — b — b Yukawa coupling. For the renormalization of the strong 
and Yukawa couplings, we employ the modified Minimal Subtraction (MS) scheme. The relevant 
renormalization constants in this work are defined as 

m b -> m b + Sm b , g s -> (1 + Sg s )g s 

Y b - Y b + SY b , 5Y 6 = n — 

m b 

b - + 

9, - {l + \8Z g )g^ (3.1) 

where g s and Y b denote the strong coupling and the h — b — b Yukawa coupling, b and g^ denote 
the fields of bottom quark and gluon. The explicit expressions of these renormalization constants are 
presented in Appendix A. 

Then the renormalized amplitude for virtual corrections, M v , can be divided into two parts and 
expressed as 

M v = M loop + m ct ^ (3.2) 
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where M loop is the amplitude from one loop-diagrams shown in Figs. 2-3 and M CT is the amplitude 
from the diagrams which contain counter-terms shown in Fig. 4. The expression of the M loop can be 
written in the form as 

M lo °P = Ui(ki)[fi^P L + f2^P R + hVly.PL + Upi^Pr + hk^P L + fek^PR + fwinhPL 

+ faPiufcPR + fak^fcPL + fiok^foPR + fulfihPL + fi2lfihPR] u j(pi) e ^(P2)TQ, (3.3) 

where Pl,r = 1 ^ 75 . The explicit expressions of form factors fi (i = 1 ~ 12) are presented in Appendix 
B, and M CT is expressed in Appendix A. 

The virtual corrections to the cross section can be written as 

& v (s,bg - h°b) = I dt 2ReJ2l(M v ^M% (3.4) 

with tmax = and t m i n = vr ^ l ~ * anc ^ again the summation with bar means the same operations as 
appeared in Ea. l)2.4|) . 

After the renormalization procedure, a v is UV- finite. Nevertheless, it still contains the soft/collinear 
IR singularities 



IR 



a s r(l — e) / Airfj, 



27rr(l - 2e) V a 



^ 6 



<» (f + f ), (3.5) 



where 



A 2 - 3 , 

^ = + 3^^^-!^^^. (3.6) 

o s — 3 s — m h 

The soft divergences can be cancelled by adding with the soft real gluon emission corrections, and 
the remaining collinear divergences are absorbed into the parton distribution functions, which will be 
discussed in the following subsections. 

2. Real gluon emission Corrections 
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The 0(a s ) corrections to bg — > due to real gluon emission (shown in Fig. 5) give the origin of IR 
singularities which cancel exactly the analogous singularities present in the 0(a s ) virtual corrections 
mentioned in above subsection. These singularities can be either of soft or collinear nature and can be 
conveniently isolated by slicing the bg —* h°b + g phase space into different regions defined by suitable 
cutoffs, a method which goes under the general name of Phase Space Slicing(PPS). In this paper, we 
have calculated the cross section for the 2 — * 3 process 

Hpi) + g( P 2) -> h°(k 3 ) + b(h) + g(h), (3.7) 
using the method named two cutoff phase space slicing method We define the invariants 

s = (P1+P2) 2 , i=(pi-k 3 ) 2 , u=(p 1 -k i ) 2 , 
s 45 = (&4 + k 5 ) 2 , i 15 = {pi - k 5 ) 2 , £ 2 5 = (P2 - h) 2 , Us = (k4 - h) 2 , (3.8) 

and describe this method briefly as follows. Firstly, by introducing an arbitrary small soft cutoff 8 S 
we separate the 2^3 phase space into two regions, according to whether the energy of the emitted 
gluon is soft, i.e. E§ < 5 s y/§/2, or hard, i.e. £5 > 5 s \/§/2. The partonic real cross section can be 
written as 

d R(b g ^ h % g ) = & s {bg ^ h o bg) + d H {hg ^ h%g) ^ (3 g) 

where a s is obtained by integrating over the soft region of the emitted gluon phase space, contains all 
the soft IR singularities. Secondly, to isolate the remaining collinear singularities from a , we further 
decompose a into a sum of hard-collinear (HC) and hard-non-collinear (HC) terms by introducing 
another cutoff 5 C named collinear cutoff 

a H {bg h°bg) = a RC {bg h°bg) + 0^(65 -» h°bg). (3.10) 

The HC regions of the phase space are those where the invariants £15, £25 ^45 become smaller in 
magnitude than 5 c s, in collinear condition, while at the same time the emitted gluon remains hard. 
<r HC contains the collinear divergences. In the soft and HC region, a s and <7 HC can be obtained by 
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performing the phase space integration in d-dimension analytically. In the HC region, a is finite and 
may be evaluated in four dimensions using standard Monte Carlo techniques [ID]. The cross sections, 
a s , <r HC and <r HC , depend on the two arbitrary parameters, 5 S and 5 C . However, in the total real gluon 
emission hadronic cross section a R , after mass factorization, the dependence on these arbitrary cutoffs 
cancels, as will be explicitly shown in Sec. IV. This constitutes an important check of our calculation. 
In the next two subsections, we will discuss in detail the soft, hard-collinear gluon emission. 

2.1 Soft gluon emission 

In the pi + P2 rest frame, the emitted gluon's (i-momentum (d = 4 — 2e) can be parameterized as 



/C5 = £5(1, ...,sm#i sin sin 0i cos 62, cos 61) 



(3.11) 



The soft region of the bg —* h b 



+ g phase space is defined by 




(3.12) 



In the soft region, the three body phase space can be factorized as[0] 




(3.13) 



with 




(3.14) 



In the soft limit, the matrix element squared for the real gluon emission, ^2\M R \ 2 , can be factorized 



into the Born matrix element squared times an eikonal factor $ e ik 




(3.15) 



where 



<Z>eik 



= N 



P1-P2 1 Pi ■ k 4 



+ N 



P2 ■ h 



(3.16) 



ipi ■ k 5 )(P2 ■ h) N (pi ■ h)(h ■ k 5 ) 



(P2 ■ h)(h ■ k 5 )' 



S 



with N = 3. The partonic differential cross section in the soft region can be written as 



da b = {Aixa s \i A r 



2e\ 



4vry r(l - e) 1 



da I dS<$> pi 



s J T(l - 2e) 2(2vr) 2 
After integration over the eikonal factors, the differential cross section is 



eik ■ 



with 



A s 



Af 



A s 



da s = da 
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a B r(l-e) / 4ttm 
2vr T(l - 2e) V s 



2\ e 



A s A s 



- — In 5 S — 3 In =■ + — In 



s — mr 



3 § — ml ' 



— In 2 5 S + 6 In 5 S In + - In 2 ■ 



— — In 8 S In 7y 

3 s — m h 6 

2.2 Hard collinear gluon emission 



s-m- h 2 
-ln 2 ^. 



s — mi 



s — m h 3 



lLi 2 [ 



s — mi 



+ 3Li 2 [- 



-u 



s — mi 



(3.17) 



(3.18) 



(3.19) 



In the limit where two of the partons are collinear, the three body phase space is greatly simplified. 
And in the same limit, the leading pole approximation of the matrix element is valid. According to 
whether the collinear singularities are initial or final state in origin, we separate <r HC into two pieces 



A-HC _ -HC , ~HC 



a 



Oa + a 



f ■ 



(3.20) 



af° is the cross section arising from the case that the emitted gluon is collinear to the initial partons, 



< ^15,^25 < S c s. And Oyf arises from the case that the emitted gluon is collinear to the final parton, 
< ^45 < S c s. We will treat <r¥ c and af^ c in Sec. 2.2.1 and Sec. 2.2.2 respectively. 



2.2.1 Collinear to the final parton 



Let and p$ be collinear to each other, < £45 < 5 c s. In the collinear limit, the three body phase 
space in d = 4 — 2e time-space dimensions may be written asjH] 

(AnY 



dT 



3 \coll 



dT 5 



16vr 2 r(l 



-dzds^[s^z{l — z)} e . 



(3.21) 



The squared matrix element can be factorized as 

y\Mf c (bg - h°bg)\ 2 ~ Y\M\bg - ^°6)| 2 e)g 2 ^^, (3.22) 

where Pbb(z, e) is the d-dimensional unregulated (2 < 1) splitting function related to the usual Altarelli- 
Parisi splitting kernels (IT]. Pbb(z,e) can be written explicitly as 

P bb (z,e) = P bb (z)+ePl >b (z), 
Pbb(z) = C F ^^, Pi b (z) = -C F (l-z), (3.23) 

with Cp = 4/3. After integration over the collinear gluon degrees of freedom, the cross section <5"j C 
can be written as [2] 



daf c = da 



a s r(l — e) f Airfj, 
2nF(l - 2e) 



1 1 A b ^ b 9 



+ A 9 , (3.24) 



where 



A[-* b9 = C F (3/2 + 2hiS s ), 

A b~>bg = c F [7/2-7r 2 /3-ln 2 8 s -lnS c {3/2 + 2ln8 s )}. (3.25) 

2.2.2 Collinear to the initial parton 

Let the hard gluon be emitted collinear to one of the incoming partons, < ti§ < <5 c si2 or 
< ^25 < &cSi2- In this region, the initial state partons i(i = b,g) is considered to split into a 
hard parton i' and a collinear gluon, i — > i'g, with p# = zpi and = (1 — -z)p«- The matrix element 
squared for bg — » /i°6(7 factorizes into the Born matrix element squared and the Altarelli-Parisi splitting 
function 

Y \M? c {bg - ^°6 3 )| 2 ~ (4™ s/ ^)V |M> 5 - h°6)| 2 f ~ 2P ^ e ) + zH^klfifA , (3 . 26) 

^ ^ \ Zt 15 Zt 2 5 J 



where 



P gg (z,e) = P gg (z)+ePgJz 



99 v 

P gg (z) = 2N[ T ^-+ l —^ + z{l-z)l P' 99 (z) = 0. (3.27) 
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Using the approximation pi — ~ zpi(i = 1,2), the three body phase space may be written as 



dr. 



3 1 coll 



dT 2 



167r 2 r(l - e 



■dzdt i5 [-(l - z)t i5 ]~ e , (* = 1,2). 



(3.28) 



Note that the two body phase space should be evaluated at a squared parton-parton energy of zs. 
Therefore, after integration over the collinear gluon degrees of freedom, we obtain^] 



da 



HC 



da 



s.0 



a s T(l — e) / 47r/i 



27iT(l- 2e) 



,2\ e 



1 



( — e)G b / A {x 1 lz)G g / B [x2) 



+ P g g(z,e)G g / A (x 1 /z)G b / B (x 2 ) + (%i ^ x 2 )\ — ( y t dx 1 dx 2 . 



(3.29) 



z z 

In order to factorize the collinear singularity into the parton distribution function, we introduce a 
scale dependent parton distribution function using the MS convention: 



1. 



Gi/A(x,Hf) = G i/A (x) + ( — ) 



a s T(l — e) / Airfi 



dz 

—Pii(z)G i/A (x/z), (i = b,g). (3.30) 

z 



2vrr(l - 2e) \ j 

By using above definition, we replace G g ^/^B i n Eq. (|3.29|) and the expression for the initial state 
collinear contribution at 0(a s ) order is 



da 



HC 



where 



da 



a s T(l — e) ( Att/j, 
27rr(l - 2e) 



,2\ e 



{G g /A{ x u^f)G h / B {x 2 ,Hf) + G g / A {xi,Hf)G b / B {x 2 ,iMf) 



+ 2J [~ h Ag c (a ^ a5)]G 9 / A (xi,^ / )G b / B (x 2 ,/i/) + (^i ^ x 2 )}dxi(ix2(3.31) 



a=g,b 



And 



with 



Af(6 - 65) 
Af (5 - 55) 



/ISC 

^0 



G a /A,B(x,Vf) 



1-Ss 



dy 



C F (2ln5 s + 3/2), 
2iVln5 s + (lliV - 2n/)/6, 
Afln(A). 



G a /A,B{x/y,Hf)P aa (y), (a = g,b), 



~ 1 — y & 

Paaiy) = P aa M$c 2) ~ P aa(yJ> (« = 



(3.32) 



(3.33) 



(3.34) 
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We can observe that the sum of the soft (Eg, ()3.18j0 . collinear ( Eg . (|3, 24|) . ()3 .3 1 f) ) . and ultraviolet renor- 
malized virtual correction (Eq.(j23J)) terms is finite, i.e., 

A% + A% = 0, 

Af + AX + A h {* h9 + A?{b^bg)+A?(g^gg)=Q. (3.35) 

The final result for the 0{a s ) correction consists of two contributions to the cross section: a two-body 
term and a three-body term <j( 3 ). 



= EL 



o r J dx 1 dx 2 da°{G g/A (x 1 ^ f )G b/B (x 2 , t i f )[Ai + A% + A h ^ h9 + A s c (b - bg) + A${g - <?<?)] 

+ G g/A (x 1 ,fi f )G b/B (x 2 ,Hf) + G g / A (xx,Hf)G b / B (x2,Hf) + (an <-> x 2 )}. (3.36) 

And 

cr (3) = J dx 1 dx2[G g / A (x 1 ,Hf)G b/B (x2,Hf) + (»i ^ x 2 )]da ( - 3 \ (3.37) 
with the hard-non-collinear partonic cross section given by 

dfr( 3 ) = -J-/ V|M 3 (6 5 -»■ /i°6 5 )| 2 dr 3 . (3.38) 
Finally, the NLO total cross section for pp(or pp) — ► bh° + X is 

CT 7VLO =CT + CJ (2) +CJ (3)_ (3>39) 

IV. NUMERICAL RESULTS AND DISCUSSION 

In the following numerical evaluation, we present the results of the cross section for the Higgs boson 
production associated with a single high-pr bottom quark via subprocess bg(bg) — ► h°b(h°b) at the 
Fermilab Tevatron and the CERN LHC. At the LHC, the &-jet is required to have a transverse 
momentum cut pr(b) > 30GeV and a rapidity cut \rj(b)\ < 2.5. At the Tevatron, the b tagging regions 
are taken to be \r](b)\ < 2 and pr{b) > 15GeV. The SM parameters are taken as: m t = 174.3 GeV, 
mz = 91.188 GeV, mw = 80.419 GeV and (Xew = 1/128 12 . For simplicity, the renormalization and 
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factorization scales are taken as /x r = nj- = mh. We use the one-loop formula for the running strong 
coupling constant a s with a s (raz) = 0.117. 

The relevant MSSM parameters in our calculation are: the parameters Mq ^ ^ and At t b in squark 
mass matrices, the higgsino mass parameter fi, the masses of the gluino nig and the A Higgs boson 
rriA, the ratio of the vacuum expectation values of the two Higgs doublets tan/3. The squark mass 
matrix is defined as 



M\ = { < ^ ) (4.1) 



a q m q m gR 



with 



m~ L = M~ +m 2 q + m|cos2/3(lf - e q sm 2 8 w ), 

rj OOO O 

rriq R = M^jj ^+m q + m z cos 2j3e q sin 6 W 
a q = A q - /i{cot/3,tan/3}, (4.2) 

for {up, down} type squarks. /| and e q are the third component of the weak isospin and the electric 
charge of the quark q. The chiral states qi and (Jr are transformed into the mass eigenstates q± and 

( Qi \ Rq ( Ql\ Rq = ( cos6 q sme q \ 
Then the mass eigenenvalues m qi and m q2 are given by 

"o 1 m? ) = R « M K R ^ ( 4 - 4 ) 
For simplicity, we assume Mq = Mq = = At = Ab = 600 GeV, \i = 200 GeV, and tua = m g = 
200 GeV by default unless otherwise stated. 

In our calculation, we use the program FeynHiggsFast |13| to generate the mass of Higgs boson(/i°) 
and the mixing angle a, and use the CTEQ5L parton distribution functions ^3]. The MS bottom quark 
mass nib can be evaluated by using the next-leading order formula JS]- In the following equations, 
we use mb(Q) to denote the MS bottom quark mass. 

m b (Q) = U 5 (Q,mb)rn b (mb), for Q < m t , 
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rn b (Q) = UQ(Q,m t )U 5 (m t ,m b )m b (m b ), for Q > m t , (4.5) 
where fn b = m b (jn b ) = 4.3 GeV. The evolution factor Uj(f = 5, 6) is 

tt in n\-( a ^) \ dU) r a 8 (Qi) -a s (Q2) tW] 

,(/) _ 12 (/) 8982 - 504/ + 40/ 2 

33-2/' 3(33 -2/) 2 1 j 

In addition, we also improve the perturbative calculations through the following replacement |15j 

— tn\ m b(Q) ( A 7 \ 



where 



2a 

Am b = -^-rrig tan (3I(m~ bi ,m~ b2 ,m~ g ), (4.8) 



with 



CL ,00 ^ 00-. C 



/(a, 6, c) = (q2 _ &2)(fe2 _ c2)(q2 _ <;2) (a 2 b 2 log ^ + fr 2 c 2 log ^ + c 2 a 2 log (4.9) 

Fig.6 shows that our NLO-QCD result does not depend on the arbitrary cutoffs 5 S and 5 C by using 
the two cutoff phase space slicing method. The two-body (c^ 2 - 1 ) and three-body (ct^ 3 )) contributions 
and the NLO cross section (a NLO ) are shown as a function of the soft cutoff 5 S with the collinear 
cutoff 5 C = 5 s /50, and tan/? = 4 . We can see the NLO cross section a NLO is independent of the 
cutoffs. In the following numerical calculations, we take 5 S = 10 -5 and 5 C = <5 s /50. 

Fig. 7 shows the dependence of the relative NLO-QCD corrections of the process pp{oi pp) — > 
bg(bg) — > h°b(h°b) + X on the mass of ^4° Higgs boson uia at the LHC. There we take tan/3 = 10, 40 
and define the relative correction as 



a NLO - a 



(4.10) 



For tan/? = 10, the relative NLO-QCD correction is about 58% to 60%. And for tan/3 = 40, the 
relative NLO-QCD correction increases from 66.5% to 73.5% with the increment of from 200 GeV 
to 800 GeV. 
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To compare the production cross section of the processes pp(or pp) — > bg(bg) — > h°b(h°b) + X in the 
MSSM with the corresponding cross section in the SM, we introduce the ratio of <jmssm and <jsm- 
The cross section ctmssm in the MSSM and the corresponding asM in the SM including one-loop 
order QCD corrections, can be expressed as 

&MSSM = <?MSSm(1 + AmSSm) 

°sm = a° SM (l + Asm)- (4.11) 

In Fig.8, we depict the dependence of the ratio (Jmssm/&sm on A Higgs boson mass at the 
Tevatron. We see that the ratio ctmssmZ&sm decreases with the increment of tha- When tua is in 
the range between 200 — 300 GeV, the cross section of the process pp — > bg(bg) — > h°b(h°b) + X in 
the MSSM at the Tevatron is greatly enhanced over that in the SM, but when vtia increases from 
300 GeV to 800 GeV, the differences of the SM and the MSSM cross section is becoming smaller. 
That is because in the MSSM parameter space for vtla = 200 — 300 GeV, the coupling factor between 
Higgs boson h° and bottom quarks sin a/ cos/3 > 1, while in the region of tua > 300 GeV, the value 
of sin a / cos (3 is getting close to 1 . 

Fig. 9 shows the cross sections of the process pp — > bg(bg) — > h°b(h°b) + X at the LHC and 
— ► bg(bg) — ► h°b(h°b) + X at the Tevatron as the functions of tan/3. We scale the mass of the 
Higgs boson rrih under the X-axis to show the changes of the phase space with tan (3. The figure shows 
that when tan/3 varies from 5 to 10, the cross sections drop quickly and m h o increases faster than 
in other tan/3 regions. But when tan/3 is larger than 10, the variation of mass of Higgs boson(/i°) is 
not sensitive to tan (3, and the cross sections drop smoothly with the increasing of tan (3. When tan (3 
approaches to the value of 50, the relative NLO-QCD corrections take the value about 70% at the LHC 
and 85% at the Tevatron. The curve behavior shown in Fig.9, comes from the feature of the phase 
space variation versus tan/3. In Fig. 10 we present the dependence of ctmssm / '<?sm on tan/3 at the LHC 
and Tevatron. We find that the cross section of the processes pp{oi pp) — > bg(bg) — > h°b(h°b) + X in 
the MSSM is greatly enhanced compared with that in the SM in all the plotted range of tan (3 = 4 ~ 50 
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with the parameter space we take. That is due to the stronger coupling strength between Higgs boson 
h° and bottom quarks with large tan f3 in the MSSM than in the SM. Then the production of the 
processes pp{ov pp) — ► h°b(h°b) + X at hadron colliders could be enhanced significantly in the MSSM. 

To study the decoupling behavior of SUSY QCD correction, we assume all SUSY mass related 
parameters have the same quantity and are pushed to a large value except \x is taken to be 200GeU. 
We denote Mq = Mjj = = A t = Ab = ijia = rrig collectively by M s and define the relative SUSY 
QCD correction as 

Sa SQCD 

&SQCD = a0 , (4.12) 

where 5a s ® CD is the cross section correction contributed by the SUSY QCD diagrams shown in Fig.3. 
In the mass region of = M s > 500 GeV, sin a/ cos/3 is close to 1. Therefore, with this parameter 
choice, the dependence of AasQCD on M s can demonstrate the decoupling behavior of SUSY QCD 
correction in the processes pp{oi pp) — ► h°b(h°b) + X at hadron colliders. In Fig. 11, we depict the 
relative SUSY QCD correction A$qcd as the functions of M s at the LHC. From the figure, we can 
see that for tan (3 = 10 the relative SUSY QCD correction is approaching to zero with M s increasing 
to 2.2 TeV, but for tan (5 = 40 the relative correction decreases to about 4% with M s increasing to 
3 TeV. That shows obviously that the large tan (3 enhances the SUSY QCD corrections and delays 
the decoupling feature. 

To analyze the scale dependence of the cross sections, we introduce the ratio of the cross section 
at scale Q to the cross section at scale Q = m h o and depict the a(Q)/a(Q = m h o) as a function of 
Q/m h o at the Tevatron in Fig. 12. The scale variation of the NLO-QCD cross section may be serves as 
an estimate of the remaining theoretical uncertainty of the high order corrections. Fig. 12 shows that 
it is evident that the one-loop NLO-QCD corrections reduce the LO scale dependence. 

In summary, we have computed the production of Higgs boson h° associated with a single high-p^ 
bottom quark via subprocess bg(bg) — ► h°b(h°b) including the NLO-QCD corrections in the MSSM 
at the CERN LHC and the Fermilab Tevatron. We find that due to the enhancement of the Yukawa 
couplings strength of the down-type quarks with Higgs bosons at large tan (3, the cross section of the 
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pp — > h°b(h b) can reach hundreds of fermi barn at the LHC and the cross section of the pp — ► h°b(h°b) 
can reach dozens of fermi barn at the Tevatron. The NLO-QCD correction reaches 50% ~ 70% at the 
LHC and 60% ~ 85% at the Tevatron in the parameter space we have chosen. 
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Appendix A 

In Appendix A, we list the explicit expressions of the renormalization constants. By using the MS 
scheme the renormalization constants defined in Ea. H3.1|) are expressed as: 

47T 

= -— GpA, 

m b ^ 

SZq = -^±(±N - -Po)A, (A.l) 

where A = \ - lE + log(4vr), /3 = (UN - 2n/)/3, N = 3 and C F = 4/3. 
We define the following symbols 

c gbb = 2 6Zg + 5Zb + 52 9 * ' 

C hbl = 5Z b + 5Y b , C bb = 5Z b . (A.2) 

Then the amplitude from the diagrams which contain counter terms M CT can be written as 
m ct = 9^.^^^^^ 

+ ^Ui{k A )[i{C gbb + ~ h) ~ C bb {i>i ~ hfhM(Pi>MT?y (A.3) 

Appendix B 



17 



In Appendix B, we present expressions of the non- vanishing form factors in Eq. (|3.3|) , We denote 
the couplings between gluino , quark and squark as 



9-ii-i'- ivgP L + vgP R )T\ g-k-b: (V^P L + vgP R )T«, (B.l) 

where 



V { } ] , = -iV2g s sin0r~ h , V { } ] , , = iV2g s cos6r~ h , 

giit.gfeift y *>&' gt 2 t,gb 2 b ys *>&' 

K (2) .j. . = -i\/2o s cos6»,-r, K (2) , = -«'\/2g s sin6>~r. (B.2) 

gtxt,gbxb ^ gt 2 t,gb 2 b a *. b v ; 

And the couplings between the Higgs boson and two quarks(squarks) are denoted as 

h-b-b: V hbb , h-k-l): V hii% ., h-ti-Vj- V hki ., (B.3) 

where the explicit expressions of the Vhbb, Vhb b- anc ^ ^hutj can ^ e f° un d i n Ref.jlbj. For simplicity, we 
introduce the following abbreviations, 

yW* . yW = pi y(l> . y<$ = p 2 

V gkb V gb 3 b V V gbib gbjb *V 



y{2)* _ T/ (l) _ ^3 T/ (2)* _ T/ (2) _ ^4 

and 



V ~ ■ V { ~> = F ■ V { ~>* ■ V { -> = F ■ (B 4) 

~gbib r V> V gbib V gb 3 b **3> 





= B o- 


B x [p^ 0,0], 


B&,Bl = Bo,B 1 [ka-Pi,0,0], 




B 0i B l 


= B o- 




0, 0] , B$, Bf = B ,B 1 [k 3 - Pl ,mg,m k ], 




B Q-, B 1 


= Bq. 


#l[-.Pl ~P2, 


m~ b ,,m § ], C§, Cg, = C , C fe /[-A;3,pi, 0,0,0], 




/~ib /~tb 


= Co, 


Ckl[-P2,-Pl 


, 0, 0, 0], C C , C c kl = C , C w [-P2, A: 3 - Pl , 0, 0, 0], 




/id rfd 


= Co, 


Cki[h, -Pi - 


■p 2 , 0,0,0], Cq,C^ = C ,C w [-A;3,pi,mt .mi. 


, m g], 


/"»/ /-i/ 
°0> U fcZ 


= Co, 


Ckl[-P2,~Pl 


,m~g,m~g,m~ bi ], C5 7 , = C , C k i[-p 2 , -pi, m~ b . 


,m~ h ,m~g], 


^/i /~ih 
^0 i 


= Co, 


Cki[-P2,h - 


pi,mg,m g ,m- w ], C™, Cj^ = C ,C k i[-p 2 ,k 3 - 


P^ m l^ m b^ m 9^ 


/~ytc fin 


= Co, 


Cki[k 3 , -pi - 


p 2 ,m- b _,m~ b .,mg], Dq,D^ = D , D k i[-p 2 ,k 3 , - 


-Pi, 0,0, 0,0], 




= D 


,D k i[-p 2 ,k 3 ,- 


- Pl ,rn~ bj ,m~ b .,rn~ b .,rn~ g ], 




Dq, Dkl 


= Do 


, D H [k 3 , -pi, ■ 


-P2,m- b .,m- b .,mg,mg] 


(B.5) 
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+ 
+ 



Since we neglect the bottom quark mass throughout the calculation except in the Yukawa couplings, 
the form factors /i, 2, 7,8,9,10 ar e irrelevant to our results and we do not present them here. 
For diagrams shown in Fig.2(SM-like QCD corrections), we can write the form factors as 

f 3 = U = [C c 22 - C c 23 + (D& - 2D a n - 2D a 21 + D a M )m 2 h + (£ft + D 25 )s 

(D a 12 -D a 13 + D^-D b 25 )t\ (B.6) 

h = k = ^fl 26 + 4Bg-32flg-16B? -120C^ + 16(2C^-Cf 2 )m 2 h 
{2C^ - 9C C - 4C a - 18Cfi + 9CT 2 + 16Cf 2 + 40C 2 C 2 - 40C 2 C 3 + 74Df 7 )t 
(18Df 2 - 29D^ + 2£>? 3 - - 6£ 2 °2 + + 2D 2 a 5 + 4D 2 a 6 )m^ 

+ (161% - 16£>? 2 + 2^ + 6D b 24 + 16Df 5 - 13D&)Sf 

+ (161% - 161% + 21% + 41% - 161% + 161% - 4Z%)f 2 ] (B.7) 

hi = /i 2 = ^^[(24 - 32B c -16B 1 c -104C 2 fe 4 + 16(2Cf 1 -^ 2 )m2)(l + s 7£) 

+ (161% - 181% - 21% + 9Z% + 91% - 18Z%)s 

+ (18Cg - 18C C - 9C?i + - 7Cf 2 - 16Cf 2 + 16C? 2 + 16C 2 C 2 - 16C 2 C 3 + 16C* 2 b 3 + 54Z%) 

+ (361% - 271% - 5Z% + 91% - 181% + 9Z%)m£ 

+ (16i% - 36Z% - 21% + 91% - 9Z% + 91% - 9Z%)£] (B.8) 

For diagrams shown in Fig.3(SUSY QCD corrections), we can write the form factors as, (The 
summation over the particle indices of squark = 1,2) has been omitted.) 

h = h( F u -> F u F fi -> F ji) 

h = 9^H [{B " ~ 8B * ~ 18C, 24-2C 2 1 + (9C& + CS- 9C 2 h 3 -C%)t)F±V hbb 
+ (8C e - (D b + 9D C + D b n + 91% + D^m-gFf^z] 



fe — fb{Fu ~^ F u Fji ~^ Fji) 
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Ai = ^^|[*(9 + 16^ e + 16^-36^ + 4Ci A + 18C[ 2 s + 18C*a - lS^mf)^^ 

+ s(9 + 165# + 16B? - 36C 2 \ + 4C l 2A + 18Cf 2 t - 18C£ 3 t - 18C fc m?)^^ 

- 2(8C e s + 8C "i - 9D c St> § 4^.- 6 J 

/12 = fuiFt^Fl F*<-i$) (B.9) 

In our paper we adopt the definitions of the one-loop integrals in Ref. ^7]. The numerical 
calculation of the vector and tensor one-loop integral functions can be traced back to four scalar loop 
integrals Aq, Bq, Co, -Do as shown in Ref.|18j. Many of the integral functions contain the soft and 
collinear IR singularities, the formulas to calculate these integrals can be found in Ref. |19j . 
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Figure Captions 

Fig.l Leading order Feynman diagrams for the subprocess of bg — ► h°b. 
Fig. 2 Virtual one-loop Feynman diagrams of the SM-like QCD corrections. 
Fig. 3 Virtual one-loop Feynman diagrams of the SUSY QCD corrections. 
Fig. 4 Feynman diagrams that contain counter-term. 
Fig. 5 Feynman diagrams for the real gluon emission. 

Fig. 6 Dependence of the cross sections for the h b production at the Tevatron on the cutoff 5 S 
with S c = 5 s /50. 

Fig. 7 The dependence of the relative NLO-QCD correction on the with tan/3 = 4, 10 at the 
LHC. 
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Fig.8 The dependence of (Jmssm/^sm on the tua with tan/3 = 4, 10 at the Tevatron. 
Fig.9 The dependence of the cross sections on tan (3 at the LHC and the Tevatron. 
Fig. 10 The dependence of ctmssm/^sm on tan/3 at the LHC and the Tevatron. 
Fig.ll The dependence of Asqcd on M s at the LHC. 

Fig.12 The variation of the a(Q)/a(Q = m h o) with the ratio Q/m h o at the Tevatron. 
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